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Abstract 



o 

li . We consider a model with non-minimal derivative coupling of in- 

1"^ ' flaton to gravity. The reheating process during rapid oscillation of the 

CO ■ inflaton is studied and the reheating temperature is obtained. Behav- 

iors of the inflaton and produced radiation in this era are discussed. 

7' 



1 Introduction 



To solve some problems of standard model of cosmology and particle physics, 
^T) ' such as the horizon problem, flatness, isotropy and homogeneity of the uni- 

^ ■ verse, the absence of magnetic monopoles and so on, and to make the bing 

bang cosmology more consistent with astrophysical data, the inflation the- 
\f^ I ory, which considers an epoch of accelerated expansion for the early universe, 

was introduced [Ij. One straightforward way to describe this era, is to con- 
cn . sider a slowly rolling scalar field , dubbed inflaton, whose energy density was 

^ I dominated by its potential during inflation [21 [3] . At the end of inflation the 

universe was cold, so there must be a procedure trough which the scalar field 
decayed to particles which became thermalized and reheated the universe. 
This could be realized by decaying of the scalar field, e.g. during coherent 
oscillation in the bottom of the potential [Ij. 
C^ ■ At first sight, it seems that a candidate for this scalar field may be 

the Higgs boson. But parameters of the standard model do not agree with 
those required for inflaton [3j. So to reconciliate the slow roll inflation with 
standard model parameters, a framework in which the Higgs boson is non 
minimally coupled to Ricci scalar has been introduced in [5j. 

Recently a model comprising a non minimal coupling between the deriva- 
tives of the Higgs boson and Einstein tensor has been proposed, which be- 
sides its capacity to explain the inflationary phase, is also safe of quan- 
tum corrections and unitary violation problem [6j. In this context, the non- 
minimal derivative coupling may allow the model to describe the acceleration 
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as well as the super-acceleration of the universe [§]• Coupling the Einstein 
tensor to kinetic term of inflaton, as we will see later, enhances the gravita- 
tional friction during slow-roll and this allows us to consider more general 
steep potentials, such as Higgs potential, without contradiction with the 
CMB observations or collider experimental bound [6]. 

In [7], this model was employed to study the natural inflation, where 
the inflaton is assumed to be a pseudo-Nambu-Goldstone boson. In this 
framework, the global shift symmetry is broken at a scale /, giving rise 
to the inflaton mass. For small field values, the potential is stable against 
radiative corrections, but slow roll requires that / becomes much larger than 
the Planck scale, giving rise to eta problem. In [7|, it was shown that non- 
minimal derivative coupling allows to take / <C Mp^ without introducing 
new degrees of freedom, and protects the tree-level shift invariance of the 
scalar field as well as the perturbative aspects of the theory. 

Similar models including nonminial derivative coupling between a scalar 
field and gravity have also been used to study the late time evolution of the 
universe by considering the scalar field as the dark energy [9]. 

In this manuscript, following [6] we assume that the inflation is imple- 
mented by a scalar field with non-minimal derivative coupling to gravity, 
with a power law potential, and study the reheating process in this model, 
which to our knowledge, was not studied before. First, we review briefly 
the inflationary epoch and then study quasi-periodic motion of the infla- 
ton at the end of the slow roll. We consider the decay of the scalar field to 
ultra-relativistic particles (radiation) via a phenomenological source during 
coherent rapid oscillation and find the reheating temperature. 

2 Non-miniinal derivative coupling model and in- 
flation 

An action describing a scalar field coupled non-minimally to gravity via its 
kinetic term is given b}U [6j: 

S = j (^ - Ig^^'d.^d,^ + ^G'^'^d.^d^ip - Vi^)) V^d\ (1) 

where G'^'^ is the Einstein tensor, ti; is a constant with the dimension of 
inverse mass squared, and Mp = 2.4 x lO^^GeV is the reduced Planck mass. 
In the absence of terms containing more than two time derivatives, additional 
degrees of freedom are not produced in this theory. 

In the spatially flat Friedmann-Lemaitre- Robertson- Walker (FLRW) space- 
time, the Friedmann equation is 



^We use units h — c — I though the paper. 



and the scalar field equation of motion is given by 

(1 + 3wH'^)ip + 3H{1 + 3wH'^ + 2wH)(p + V'{ip) = 0, (3) 

where "dot" denotes derivative with respect to time and "prime" denotes 
derivative with respect to ip. 

Using the energy momentum tensor derived from ([1]) , the energy density 
and the pressure of the scalar field are obtained as 

. 2 

P^ = ^-Viip)-^i3H^ + 2H)^^-2wH^yi, (4) 

respectively, p^ and P^ satisfy the continuity equation 

p^ + 3H{P^ + p^) = 0. (5) 

In the presence of another component, with the energy density p-ji and 
the pressure P-ji, interacting with the scalar field via the source term Q, the 
continuity equation becomes 

p^ + 3H{P^ + p^) = -Q, 

pn + 3H{Pn + pn) = Q- (6) 

In the inflationary era, we take the field ip as the inflaton, and assume that 
the universe is dominated by only this scalar field. But in the subsequent 
epochs, one must also take into account the presence of other components 
such as radiation. 

In the following we adopt the high friction condition [lOj 

wH^ > 1. (7) 

The slow roll regime is characterized by 

(p < 3H<p, ^ ^ ^' 

9wH^ip^ < 2V{ip), (8) 

V{ip) 



yielding 



H^ ~ 



3M|, 



V'{^) 



^ - -9^ ^'^ 



The slow-roll conditions are satisfied provided that 

V"{ip) 3w 

In comparison with minimal models, conditions (jlOp may be satisfied by 
more steep potentials in the high friction limit ([Tj) [lOj . 
For the potential 

y((^) = \^\ (11) 

these relations require ip'^^'^ ^> -^. In the slow-roll era, where ti <^ H^ 
holds, the requiring that the model be outside of the quantum gravity regime 

implies that R ~ 12H^ <^ —^. This condition, in terms of the potential, is 
rewritten as 

FM « -^. (12) 



The number of e-folds is given by 



W /-^O y2(^) 



where ipe ((/'o) is the value of the scalar field at the end (the beginning) of 
the slow-roll inflation. 

By assuming ipe <^ ipo, one can estimate the e-folds number for a chaotic 
inflation with potential ()lip as 

When the conditions (jlUp cease to be valid 

V"{ip) ^ 3w 



p 



/2 



vHp)^m. 



(15) 



the slow-roll inflation ends. For power law potential pip and q ~ 0(1), (|15p 
gives 



3 Quasi-periodic evolution 

In this part we try to study the dynamics of the inflaton after the end of 
slow-roll. In this era, conditions (jSllOp cease to be valid and quasi-periodic 
evolution of the scalar field about the bottom of the potential begins. In 
fig-©! using numerical methods, oscillation of the field ip is depicted for 
quadaratic potential 

V{^) = \m\\ (17) 

This figure shows that, at the first stages after the slow-roll, the amplitude of 




Figure 1: ^ 

with initial conditions {0(1 



■^ in terms of dimensionless time r 



0.056, 6(V\ 



mt, for {wTTT? = 10^ 
0}, for the quadratic potential. 



the scalar field drops down rapidly. Later, during a phase of rapid oscillation 
of the field, the rate of decrease of the amplitude becomes much less than the 
oscillation frequency. To get an an insight about the solutions in this epoch, 
we proceed in the same way as [11] and consider the power law potential 
(jlip . with an even integer q. 

During the quasi-periodic evolution where the amplitude and the fre- 
quency of oscillation are time dependent, the scalar field is presented as 

m 



ip{t) = $(t)cos / A{t)dt 



where, the amplitude $(t) is given by 



(18) 



(19) 



This equation means that the potential, when is evaluated at the amplitude, 
gives all the energy. A(t) is some function of time which may be determined 
as follows: By taking a time derivative of (jlSp . we easily obtain 



sm 



A{t)dt 






(20) 



20D, and ([IHD result in 



A' 



viMl 






$2 _ (^2 

The continuity equation and eqs. dH) lead to: 

p^ = -'iH[ij?{l + SwH"^ - wH) - 2wHip(p\. 



(21) 



(22) 



By substituting (p from the equation of motion ([3]) into the above equation 
we arrive at 

p^ = -3H[{9wH^ + 3wH)ip^] + 2ipV^. (23) 

By using $ = ^^ derived from (fT9|) . we deduce 






gP(^ V 2 p^ 



(24) 



where +(— ) corresponds to 99 > 0(< 0). This expression is much less than 



unity, 






^ 1, whenever the slow-roll is ceased and 



$<( 



qMJ, , 



19+2. 



(25) 



becomes valid. This is opposite to slow-roll conditions (see eqs. ()8|10p and 
their subsequent discussion). In this regime 



A' 



^' 



$2 _ ^2 



2ip^-Viip)) 

9u;i?2($2 _ y,2) 

2M|, {p^ - Vjif)) 
3wp^{^^ - (^2) 



(26) 



I « 

when 



is the stage of rapid oscillation or high frequency regime, i.e. 
<C A. In this epoch we have 



6 


= 


2H 
~qH 


= 





^ = -±L = ±^ «^, 



(27) 



showing that the amplitude, the Hubble parameter and the energy density 
decrease slowly during one period of oscillation. 

The time average of energy density over an oscillation cycle (from t to 
t + T) is 



(P^(i)> := 



1,'+^ pAt')dt' 



(28) 



where T is the period. As the ampUtude decreases very slowly, we deduce 
{Pip{t)) = y(<l>(t)). This result is in agreement with {[H]), because from ([27|) . 
it is clear that p^ changes insignificantly during a period. 

During the rapid oscillation of the scalar field, the parameter 7, defined 

by 7 = 1 + 7-^, is determined as follows 

/q tr2 -2 diwHifi'^) 

^^- — w> — '- '^'' 

To obtain the above equation we have used ^. Using ([!]), (j29|) can be 
rewritten as 

2 {p^ - Vjip)) 

3 (p^) 

2 J% V/Oy - Vi^)d^ 



2 /J VI - x'idx _ 2q 



1 1 , -o„,.- (30) 



/oTfe^- ^'^ + 6 



In the above computations, we have converted time integration to (/? inte^ 



gration, and the variable change x = ^ was applied. The same method for 



i/; = 0, results in 7 = -^ [TT]. This constant is the effective value of 7 in 
rapid oscillation era, i.e., in this era the scalar field behaves as a barotropic 
fluid. To elucidate this point we proceed as [Hj. Taking time average of ([S]), 
and by using ([30]) . we obtain 



2q_ 
q + 2' 



{p^) + zr^H {p^) = 0. (31) 



Using the definition of time average, we obtain (p^) = -^, where 6p^ is 
the change of p^p over the period T. Hence 

+ ^H{pp)=0. (32) 



5/v _2g^ 
T +9 + 2''^'^^' 



In the high frequency regime -^ w p'^, and we arrive at 



P^ + -^Hpp = 0. (33) 



In the following, as in /133\) . we will not use the symbol <> for the oscil- 
lating scalar field, e.g., by p^ we mean the time averaged of the scalar field 
energy density in the sense explained above. 



Prom eq. (|19p we get an approximate equation for <I> evolution 

6 ^ -^^H^. (34) 

q + 2 ^ ' 

From ()34p and ()27p it is obvious that H <^ A. This means that the expansion 
rate is much less than the oscillation frequency. Effectively the scale factor 

is a{t) oc t 9 and the Hubble parameter is given hy H = 2— 1~^. p satisfies 

%^ = 0, (35) 

at 

corresponding to p^ ex t~^ which is q independent. 

4 Particle production 

In this part we try to study infiaton decay to ultra-relativistic particles 
(radiation) in the rapid oscillatory phase. This decay is due to the inter- 
action between the infiaton and produced particles. To study this decay, a 
phenomenological interaction (see(l6|)), 

Q = ^vb^ (36) 

where T is a positive constant, was proposed by [T2]. Afterwards, in [13) . 
a Lagrangian, including bosonic and fermionic fields and their interactions 
with the infiaton, was introduced, and it was shown that the effect of par- 
ticle production can be explained by adding a polarization operator to the 
infiaton mass term. There was shown that, for quadratic potential, the role 
of polarization operator may be mimicked by the phenomenological friction 
term ()36p . in high frequency regime. However, by taking under consideration 
the back reaction of the quantum effects on the evolution of the infiaton field 
[14j . and also considering possible decays of the infiaton to other particles 
[IS], ([36]) can no more be deduced from a Lagrangian. 

Other phenomenological models with temperature and field dependent 
friction coefficient term, Q = T{T,ip)ip'^, have been also considered in the 
literature (16j . 

However, as the nature of the infiaton and also primordial produced 
particles are unknown, we have not yet an exact expression for the form 
of interaction. In our study, we adopt the widely used phenomenological 
interaction ()36p . which reduces significantly computational complexity. 

In the presence of the source term ()36p . the infiaton evolution is given 

by 

3wH^^ + 3H{3wH^ + 2wH)ip + V'{(f) = -Ttf. (37) 



In the following we consider the power law potential (jlip . Using (j3]) and 
with the same method used in (j30p . one can obtain 



{9wH^^'') = -^P^. (38) 

Inserting the above relation in the first equation of ([6]), and by replacing 
Pifi + Pip with its time average over an oscillation cycle, we arrive at 

7r 

p^ + 3H-fp^ + -^^^Pip = 0. (39) 

We note again that in the above equation, all the values must be regarded 
as the time averaged values over one oscillation. This relation is only valid 
on large time with respect to the period of fast oscillation. The second 
term (friction term) describes dilution resulted from the universe expansion, 
while the third term corresponds to particle production during the coherent 
oscillation of the inflaton. 

Comparing ()39p with the corresponding minimal model relation 

p^ + 3H-frnPifi + Im'^Pv = 0, (40) 

where 7^ = t^, shows that, in the limit wH"^ » 1, the decrease of pip 
(due to particle production) and as we will see the reheating temperature 
are very less than the corresponding values in the minimal model. Note ()39p 
is true only for wH^ ^ 1, so (|4Up cannot be derived from (|39p by simply 
setting u; = 0. 

The solution of (l39l) is 



Pip (X a '^ exp 



ill 

^W -I to. 



t H-\t')dt' 



(41) 



where tosc. is the time when the oscillation commences. In the beginning of 
particle production, the universe is dominated by the scalar field, and we 
assume T <^ wH^ , but H is decreasing and this approximation ceases to be 
valid later, when the second and the third term in (|39p acquire the same order 
of magnitude (we will denote this time by trh)- So, with our assumptions, 
and in the scalar field dominated era it is safe to use the approximation 
H = Yj until t = trh- The behavior of the Hubble parameter, during rapid 
oscillation, can also be investigated via numerical method. E.g. in fig. ([2]), 



H is plotted for the quadratic potential (|T7|l and by using Friedmann and 
continuity equations. So let us write (jH]) as 

p^ = Et-'exp{-'^t^), (42) 

where S = t^^^ /9osc.e^™'*°'"=- and posc. = p{tosc.)- The decrease of p^p is due to 
particle production which is encoded in the exponential term and also due 
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Figure 2: h := — in terms of dimensionless time r = mt, for \wm? = 
IqS, ^ = 10-2,}, with initial conditions {(/>(!) = 0.056, (/>(1) = 0}, for the 
quadratic potential. 

to the term t^^ corresponding to dilution via the universe expansion. The 
dilution term is independent of q as was explained after eq. ()35p . For larger 



values of 



wM'^ 



the decay rate is faster. 



Comparing this result with what was obtained before for {w = 0, 7 
0.5} model [I7|, 



P<P 



;-y^)^exp[-r(t- tosc.)], 



(43) 



shows that the decrease rate of p^, due to the expansion of the universe, 
is the same. But in the presence of tu, the rate of particle production is 
decreased. We note again that our result ()42p is true only for wH^ ^ 1, so 
we cannot obtain (|13]) from (|12]) by setting w = Q. 
The radiation satisfies 



pn + 4:HpTi = r7 



Pv 



SwH"^' 



(44) 



To study the evolution of p-ji during scalar field dominated epoch, i.e., when 



H' 



3Mp 



p^ we write (|44p in the form 



p-R, + AHpTi = 

whose approximate solution is given by 

3r72M|, 



FtmI 

w 



8+37 



Pn 



(8 + 37)^; 



t[l_(:£££:)^^]. 



(45) 



(46) 



To derive the above equation we have assumed that after the slow-roll, 
the universe was cold: pnit = tosc.) = 0. In fig. ([3]), pii is plotted for 
the quadratic potential ()17p . showing that the radiation density increases 
monotonically in the rapid oscillation phase. However this behavior is valid 
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Figure 3: dimensionless energy density f j°^ (points) and §?,a (line) in 



1-2 



terms of dimensionless time r = mt, for \wm = 10 , — = 10 ,| with 
Pose. = 8.3 X 10~'^m^Mp, for the quadratic potential 



only for p^ > p-ji, and in the radiation dominated era we expect that pji 
decreases. 

Based on Friedmann and continuity equations, the behavior of p-ji is also 
depicted in fig.Q, via numerical method for the quadratic potential (fTTl) . 
This figure shows that the increase of p-n continues until p-ji ~ p^. After 
t = trh, i-e. from the beginning of radiation dominated era, p-ji begins to 
decrease. 
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Figure 4: dimensionless energy density 



Pv 
ImH 



p-n 



(initially upper graph) and 
' 10^ r ^ 10-2 I 



-n-itTT in terms of dimensionless time r = mt, for \wm 

with initial conditions {0(r = 1) = 0.056, (f){T = 1) = 0}, for the quadratic 
potential. 

Relativistic particles interact quickly (with respect to the expansion rate 
of the universe) with each other to become in a thermal equilibrium char- 
acterized by temperature T^, given by 



Pn 



30^*^' 



(47) 
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where 5* is the total number of effectively massless degrees of freedom. If 
Tj. is greater than the electroweak scale Tj. > SOOGey then g > 106.75 [18j . 
The reheating time, trh, is given by 

r ~ dwH^itrh), (48) 

i.e. , when the second and third terms of (|39|) acquire a same order of mag- 
nitude, as it was mentioned before. At this order of time, pnitrh) has also 
the same order of magnitude as Pip{trh)- At the reheating time we have 

H ~ (9:^)^) so the use high friction condition ([7]) during the reheating era 
is safe only when 

wT^ > 1. (49) 

For t > trh , almost all the energy of the inflaton is transferred to newly 
produced particles and the universe becomes radiation dominated. 

Now, we can estimate the reheating temperature defined by Trh = T{trh)- 
From (gU]), and ^f^ <C 1, we obtain 



rh 



Trh ^ 1.89 ( -^y g,-^Ml,{-)l (50) 

Note that the above equation could also be obtained by using H^ ~ -ttttPti, 

and ()47p . This temperature is specified only by parameters of the system and 
is independent of initial conditions. In the absence of non-minimal derivative 
coupling, i.e. for ti; = 0, and for 7 = 0.5, the radiation energy density is 
obtained as [T7] 

^^ = To^[^-(— )^]- (^^) 

Pti increases rapidly from p^i = to its maximum value and then decreases 
again, so the maximum temperature is occurred in the beginning of (/? os- 
cillation before reheating (see fig.©, depicted for the potential (fT7|) ). This 
is in contrast to our model, where as it can be seen from (06]), p-ji increases 
continuously in tp oscillation epoch until p-;^ ~ /Oi^. In u; = model, reheating 
temperature is determined by [T7] 

T,;,(u7 = 0) ~ 1.2g;^M|ri (52) 

Therefore Trh{wH'^ » 1) <C Trh{w = 0), provided that a same F is taken 
into account for both theories. Note that as the reheating process is realized 
after inflation, Trh must be below the GUT scale: Trh < W^^GeV. 

Based on astrophysical data, we are able to estimate the relation of 
reheating temperature and the number of e-folds. To do so, we proceed as 
follows: Consider a length scale I which at time t = t^,, m the inflation era, 
left the Hubble radius: 

/ = ^-1- (53) 
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Figure 5: — tttt, in terms of dimensionless time r = mt, for w 
with Pose. = 8.3 X 10~*?7i^ Mp, for the quadratic potential. 
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The number of e-folds from t.^ to the end of the sfow-roU inflation, denoted 
by subscript e, may be expressed as 



M,. 



ln(^ 



In 



0-rh 0-eq ttQ a^H^: Hq 



(54) 



where 0, rh, and eq subscripts denote present, reheating, and matter-radiation 
equality densities epochs respectively. As we have seen, during rapid oscil- 
lation (from te until trft)) the universe is dominated by a scalar field whose 
the effective equation of state parameter is given by 7 — 1, whence [19] 



K = 62 + ln (^) +ln 



v) 



lO^^GeV 



+ ilnf^Vf---V-f- 



(55) 

The right hand side expressions are determined as follows: By setting oq = 1, 
and adopting the WMAP (pivot) scale I = ^ = 500Mpc [20], we arrive at 



In 



Li* -LI J; 



In 



Ho 

k 



In 



0.7 
~6~ 



(56) 



The present Hubble parameter is Hq = 30000 Mpc ^ [20j. To obtain T4, we 
consider the scalar perturbations. The spectral index Ug is [10] 



n.o 



1 



Mj, 



w 



F2 



2V"{^) AV'{^) 
3 V{ip) 2,V{^Y 



(57) 



The quantities in the right hand side must computed at the time of hori- 
zon crossing Csk = aH during slow-roll inflation, where k is the comoving 
wavenumber and Cg is the sound speed of scalar perturbation. In the limit 
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wH^ S> 1 we have c^ ~ 1 [21], and for the potential (fTT]l . we deduce 



K = aA (H^#±^V'\ (58) 

V W{1 - Us) J 



Based on WMAP data [20], 

Ws = 0.968 ±0.012. (59) 

From t^K until ig; the slow-roll approximation is valid, hence Ve may be esti- 
mated with the help of (fTSj) . as: 

K = A^ ^i^ j . (60) 

Prh is specified by (P6|) or equivalently by (|Tr|) and ([50]) : 

'^ 8 + 37 Vw^y 

1 j. j. 

In the slow-roll regime, as expected, we have F/ ~ Ve^ . VJ' < IQ^^GeV also 

holds [19] (lO^^GeF is the GUT scale). The reheating is provided by the 

inflaton energy, hence p^h ^ K- So finally we expect to have A/"* < 60. Only 

for a prompt reheating, A^ ~ 60 may be possible. 

As an example, for the quadratic potential V{^p) = ■^'^'p'^ where 






\/2mM|, 
a/(1 -ns)w 



w 



2 



Prh = 0.15M|,(-) . (62) 



A/"* depends on the parameters of the model, i.e. F, w and m, e.g. choosing 
wm? ~ 10^, ^ ~ 10~^ (in agreement with (jl9|) ). and setting m = 10~®Mp 
[2], we find A/; = 43.41. 

5 Conclusion 

As a summary, we briefly discussed inflation in the framework of a non- 
minimal derivative coupling model proposed in [6|, and then studied a gap 
in the literature: the reheating process in this framework. We investigated 
inflaton evolution in quasi-periodic oscillation at the end of slow-roll. We 
allowed the scalar field to decay to ultra-relativistic particles (radiation) via 
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a phenomenological source term. We obtained the reheating temperature 
which was independent of initial conditions. 

We showed that the energy density of radiation, during oscillatory era 
and when it is smaller than the energy density of inflaton, increases mono- 
tonically. This behavior is in contrast to ordinary inflation theory where the 
maximum temperature occurs before reheating. We confirmed our results 
via numerical methods. 
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